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SUMMARY 


In this report the method of small oscillations is used to investi- 
gate the stability of a stratified plane laminar flow (x and y are 
rectangular coordinates in the plane; x direction = primary flow 
direction, y direction perpendicular upward) . The laminar flow U = TJ(y) 
serving as basis is a "boundary-layer flow" which increases from the 
value U = 0 at the wall to a constant value 11 = 1^ at great wall 
distances. Stratification of densities is assumed in the boundary 
layer (thickness =5), while outside of the boundary layer, in the zone 
of constant velocity TJ m , the density is assumed to be constant. The 
flow function of a partial oscillation of the superposed periodic 
disturbance motion is written as 

(x,y ,t) = cp(y)e i ( CU£: “P t ^ = cp(y )e la ^ ;x: " c ’ fc ^ 


where a is always real and identifies the spatial angular frequency of 
the disturbance (X. = 2n/a = wave length of the disturbance) ; P is, in 
general, complex, and the prefix of the imaginary part of P decides 
the stability or instability of the disturbance. The stability investi- 
gation is a characteristic value problem of the differential equation of 
the disturbance 


(U - c) 2 [cp" - a 2 cp -(U - c)u"cp + Kcp + L(U - c)jj-(U - c)<p' + U'cpj 


= -~(U - c)[q>"" - 2a 2 cp" + a^cp - L(cp" ' - a 2 9 ')j 


for the disturbance amplitude cp, with the boundary conditions cp = cp ' = 0 
at the wall y = 0, and y = » at great distance from the wall; 


*"Turbulenz bei Warmeschichtung. " Zeitschrift fur Angewandte 
Mathematik und Mechanik, Band 15* Heft 6, December 1935* PP* 313-338. 
Lecture presented at the fourth international Mechanics Conference held 
at Cambridge, England, from July 3 to July 9* 193^* 
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E denotes the Reynolds number, K = 


b dp g& 

p ^ Um 2 


- (Richardson number) 


and L = 


B dp 
P <3y 


are nondimen&ional stratification numbers 



Is a Froude number. For a specified group of four values a, R, K, L, 
this equation has only one solution for a particular value of c. The 
study is limited to pure real values of c, that is, to the determination 
of the curve of the transition points from stability to instability 
(indifference curve) for specified values of K and L in the aE 
plane. Stable (K > 0, L > 0) and unstable stratifications (K < 0, L < 0) 
are investigated. It 1 b found that for constant Froude number with 
increasing stratification (increasing K) the critical Eeynolds number is 
greater and the region of unstable disturbances in the aE plane smaller 
until finally at a critical value of the stratification quantity K, 
which still depends on the Froude number, complete stability of flow 
prevails for all disturbance wave lengths and Reynolds numbers. 


The calculation is carried out for the flow past a flat plate, for 
which To ll mi en had computed the indifference curve for homogeneous flow 

a dp 1/ dU\ 2 

The critical stratification quantity a — > 


■ ) , at which the 


p 

turbulence must become zero, ranges from 0- 0409 to 0-029 for Froude 
*m 2 

numbers = 0 to 5; while Taylor and Goldstein arrived at cL- =1/4 by 

5g 


a similar calculation in which the fluid friction and the curvature of 
the profile had been neglected, and Richardson and Rrandtl had 
obtained = 1 and ©k = 2, respectively, by rough estimates. 


The essential premises of the present investigation (plate profile, 
stratification of densities in boundary layer only) are fairly well 
confirmed by Eeichardt's measurements in the Gottingen hot-cold air 
tunnel, where the upper plate is heated with steam, and the lower plate 
cooled with tap water (stable case). The decision, whether a measured 
velocity profile was laminar or turbulent, was made from oscillographic 
records of the voltage fluctuations of a hot wire. The comparison of 
the measurements with the present theory indicates a very satisfactory 
agreement . 


1. INTRODUCTION 


On cool summer evenings, when there is a slight breeze, a person 
can occasionally observe the floating of damp fog with a very sharply 
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defined "boundary over a wet meadow, a sign of the fact that the turbulence 
of the wind has stopped completely, so that the air layers slide laminar 
over one another without turbulent intermingling- This is due to the 
development of a marked temperature difference as a result of the evening 
cooling, which prevents the warmer and hence specifically lighter upper 
layer frcm mixing with the cooler, heavier air layers near the ground. 

Closely related with it are Richardson's observations (reference 1) 
of the Influence of the gustiness of the wind which can be regarded as 
a measure for the turbulence strength on the vertical temperature 
difference. He observed that the gustiness is small when the air is 
colder below than above and that it increases when the temperature 
difference approaches the unstable adiabatic equilibrium- But, 
surprisingly, the gustiness does not Increase much more when the unstable 
superadiabatic temperature difference prevails. 

Taylor had made similar observations in 1916 (reference 2) . Rrcam 
simultaneous records of temperature and wind velocity at about 40 meters 
height above ground for several days and nights, he found that on the 
nights when the nocturnal temperature minimum was very marked as a 
result of strong radiation, the turbulence was completely gone- But on 
nights with lesser cooling, such as result from cloudiness, the wind 
velocity fluctuations were almost as great as during daytime. 

The flow of fresh water over salt water without substantial inter- 
mingling, as observed in the Kattegat, for instance, belongs to the same 
group of phenomena. Even the surprising stability of Bjerknes's polar 
fronts, where the cold air masses form a wedge under the warm masses, 
is traceable to the stratification (reference 3) • Taylor, in 1927, made 
a simple experiment by which the stabilizing effect of the stratification 
of densities on the turbulence can be shown by means of a salt solution 
(reference 4) . 

Conversely, there is Increased turbulence, and hence a stronger 
intermingling, as a result of convection motions^ and, when following 
strong radiation, the lower air layers are heated more than the upper 
layers (reference 5) • 

Reichardt has conducted experimental investigations on a stratified 
flow since 19 2 7, under the direction of Professor Prandtl. A stream of 
air in a horizontal channel Is blown between a plate cooled with water 
and heated with steam, whereby the upper, as well as the lower, plate 
can be heated (stable or unstable stratification). In 1927, when these 
experiments were still in the Initial stage, Rrandtl set up a simple 
theorem, in form of an energy consideration (reference 6) , which Is 
briefly reviewed. 
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Assume a flow in horizontal direction, with a stratification of 
densities in perpendicular direction (y-direction) , so that the density- 
decreases continuously upward. During the turbulent mixing motions, 
work is then performed by the fact that heavier matter is raised and 
lighter matter lowered with respect to the lift. The path traversed by 
a particle in vertical direction, before it mixes again with the new 
surrounding, is the Prandtl mixing path l (reference 7) , and the 
difference in lift per unit volume of a particle shifted in the vertical 


over the path length l 


is 



The work per unit volume of displaced 


fluid mass while traveling over the path length l is accordingly 


r y= ^i +2 / \ 

/ (' g<y ■ yi) i) ^ = ■# * 2 1 

J y^i v 

To identify the quantities participating on the exchange, visualize 
a level surface F} on one fraction j3j_ of this area, an upward motion 
with velocity vp prevails, and on a fraction P 2 a downward motion 
with velocity V 2 , so that the total flow volume in unit time is 
F(PiVi + 32 T 2) ’ With it, the lifting work due to the stratification of 

densities is 


» -F (^v-l + P 2 v 2 ) ~ S l 


2 

<47 


(1) 


This work must be supplied from the stored energy of the turbulent 
mixing motions. This is obtained by the work of the basic flow at the 
element, which is given by the product of the apparent turbulent shearing 
stress with the displacement velocity. For a body of base area F and 
height 2, in which the afore-mentioned volume is exchanged, this work 


is 


Ft2 


du. 

dy J 


t is the turbulent shearing stress for which t 


pu'v ', 


according to Reynolds, where u' and v' are the turbulent fluctuation 
velocities (u = u(y) = mean flow velocity). According to Prandtl, 


u' = 2 — and, according to the foregoing, v' = p^v^ + P 2 T 2' Wp a 

numerical factor K , w hich is still to be determined, in which the 
correlation factor of u'v 1 enters, the turbulent shearing stress is 
then 

t = kP^( P 1 t 1 + |3 2 v 2) 


■■11 11 


■ ll■llll■l 11 ii iiiiii 
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hence the work L-^ of the turbulent apparent friction 


L-fc = icFpl 2 




( 2 ) 


The existence and disappearance of the turbulence as a result of 
the stable stratification depends upon which of the two energy amounts 
is the greater. If the work of the turbulent apparent friction is 
greater ^L-fc > L a ^, the difference for maintaining the turbulence remains, 
but, if the work against the weight difference is greater ^L a > L-fc), the 

turbulence must die out. Turbulence is therefore possible on an energy 
basis when 


FkP ^ 2 



+ 



> 



p2 T 2)lb 2 if 


or, after abbreviation 



2 dy 


For the nondimens ional quantity 


£&P 

pdy 



which serves as basis for all flows with stratification of densities, 
the "Richardson number" is introduced, since he was the first to study 
stratified flows in 1920 (reference 9) • Thus 0=0 denotes the 
homogeneous fluid. By Prandtl's energy appraisal the result is as 
follows: 

Turbulence, on an energy basis, is possible for 0 < 2k 
Turbulence, on an energy basis, is impossible for 0 > 2k 
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Although, nothing definite can he said about the numerical factor k , it 
is suspected that it lies near unity. Prandtl chose k = 1 and so 
obtained 0=2 as stability limit. • 

Taylor indicated subsequently (reference 10) that a factor 2 cancels 
out when Prandtl's considerations are refined, so that Prandtl 's 
stability limit 0=1 would agree with similar considerations by 
Richardson (reference 9) • Taylor (reference 11) and Goldstein (refer- 
ence 12) later continued the theoretical study of flow with stratification 
of densities, that is, the stability of a plane laminar flow with stable 
stratification of densities, by the method of small vibrations- Viscosity 
and compressibility were disregarded and the velocity profiles were 
limited to those consisting of straight pieces, for reasons of mathe- 
matical simplicity. In the differential equation of the disturbance, 
only the effects of the stratification of densities on the potential 
energy (gravity effect) were taken into account, while the inertia effect 
of the stratification of densities is neglected. /The cases in which this 
is permitted, and the resulting simplification in the calculation, will 
be discussed later. The results of Taylor's and Goldstein's stability 
investigations can be represented in the two nondimensionals 0 
and l/&, 1 = wave length of disturbance, 5 a characteristic length 

of the velocity profile (boundary-layer thickness) . For the majority 
of cases explored by Taylor and Goldstein, no actual stability limit 
resulted for ©, but there still remained a certain although very narrow 
range of unstable disturbance wave lengths X for every value of ©. 

Taylor obtained a definite result for the case of a fluid extended 
infinitely upward or downward, with laminar velocity distribution and 
uniform density distribution; he obtained 0 = iA as stability limit. 
Goldstein arrived at the same stability limit for the case of a fluid 
extending to infinity upward or downward; with -uniform density distri- 
bution, the velocity below and above is constant, but varies linearly 
with the height in an intermediate layer. 

A comparison of this theoretical result with test data from the 
Gottingen warm-cold air tunnel produced no agreement at all relative 
to the stability limit 0 = l/k, where, of course, it should be borne in 
mind that the velocity profiles measured in the tunnel are not linear, 
as postulated in the Taylor-Goldstein theory. It therefore seems 
appropriate to study the modification of the Taylor-Goldstein calculations 
with due allowance for the friction and to choose for the velocity distri- 
bution of the laminar flow a profile better adapted to the conditions to 
be realized by experiment. This appeared to be all the more promising, 
as the stability study on the homogeneous flow itself produced a satis- 
factory result only when the friction was included in suitable manner. 

Tollmien (reference 13) indicated that profiles with finite curva- 
ture other than zero must be used as basis, since the approximation of a 
laminar profile by pieces of straight lines is insufficient. 
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In the first stability investigations of the laminar flow of a 
homogeneous fluid by Lord Rayleigh (reference lh) , both the friction 
and the profile curvature were disregarded; they did not give the looked- 
for instability- Further investigations by Sammerfeld (reference 15) , 

R. v. Mises (reference 16) , and L. Hopf (reference 17) on the Couette 
flow (linear velocity distribution) , which partly included the friction, 
also failed to give the desired results; it resulted in stability for 
all disturbance wave lengths and all Reynolds numbers- Subsequent 
investigations by Prandtl (reference 18) and Tietjens (reference 19) , who 
allowed for the m aximu m friction terms and used profiles consisting 
of- pieces of straight lines as basis, obtained an instability for the 
first time, but still no stability l imi t - Complete success, that is, the 
correct theoretical calculation of the stability limit (critical Reynolds 
number) was bestowed on these investigations only after Tollmien took 
the curvature of the velocity profile also into consideration. 

The Taylor-Goldstein stability studies of the flow of stratified 
fluids form the analogy with Rayleigh's study on homogeneous fluid, as 
they ignore friction, as well as profile curvature. They are to be 
extended in the following by taking friction and profile curvature 
into account. (Compare the subsequent outline.) 


Outline of Past Stability Investigations 



Homogeneous fluid 

Inhomogeneous fluid 
with stratification 
of densities 


Without 

friction 

With 

friction 

Without 

friction 

With 

friction 

Linear profile 

Rayleigh 

Sommerfeld, 
v. Mises, Hopf, 

Prandtl, Tietjens et al 

G- I. Taylor 
Goldstein 


Curved profile 


Tollmien 


Schlichting 


Closely allied with these flows with stratification of densities 
are the curved flows of homogeneous fluid, as, for example, the stable 
stratification due to the centrifugal forces in the two-dimensional flow 
between two concentric, rotating cylinders, of which the one on the inside 
is stationary, while the one on the outside rotates. In 1927 Prandtl 
obtained, by an energy study similar to that on flow with stratification 
of. densities, a stability limit for a flow stratified by centrifugal 
forces, which is in close agreement with Wendt's measurements (reference 20) 
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The effect of such a stable stratification of centrifugal forces on the 
critical Reynolds number was studied in an earlier report (reference 21 ) , 
and the results were also in satisfactory accord with experiment. 

Returning to the flow with stratification of densities, the stability 
investigation can very generally be formulated as follows: 

For each superposed disturbance (wave length h) , the magnitude of 
amplification (damping) is to be computed for each Reynolds number and 
for each layer (Richardson number), the laminar flow U = U(y) being 
specified. If A 0 is the amplitude of disturbance at time t = 0, 

P ± t 

the amplitude A at time t is A = A 0 e , where pj is a 

measure of the amplification. This problem is dependent on eight 
essential factors, namely. 


Pi 

X 

5 

U m 

P 

Ap 

v 

g 


log de ( in) cr ement of the amplification 
wave length of disturbance 
characteristic length of laminar flow 
maximum velocity of laminar flow 
density 

density difference in laminar flow 
kinematic viscosity 
gravitational acceleration 


From these eight quantities, five independent nondimens ionals can 
be formed which are the characteristic variables of the stability 
problem, namely, 

P-jS 

= dimensionless amplification quantity 


— = dimensionless disturbance wave length 
& 


u m s 


= R = Reynolds number 
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gAp 
S p 



= 0 = Bichardson number 


^ = corresponding density variation 
P 



Um _ Ap 
5 g ~ P I 


0 = Froude number 


A. 


/ 


The most general solution of the stability problem therefore 
implies the finding of the function 


G 




0 


(3) 


In an earlier report (reference 22) , this very general function G 
of five variables had been specialized to the effect that the amplifi- 
cation for a homogeneous fluid (0=0, Ap/p = 0) was computed as 
function of the Eeynolds number and the wave length of the disturbance, 
hence 


Gl 




= 0 


The subsequent investigations are based upon the limitation 
employed in almost all other stability studies, that only wave lengths 
of disturbance X = X Q which are neither damped nor amplified are 
involved (neutral disturbances) . 

Hence, it is assumed that Pi = 0 and the more special function 



(4) 


which depends on only four variables is determined in place of (3). In 
contrast with the homogeneous fluid there results at once through the 
introduction of the density stratification two new nondimensionals, 
namely 0 and Ap/p, where 0 gives the effect on the potential energy 
of gravity of the density stratification and Ap/p represents the effect 
of the stratification of densities on the inertia. The ratio of these two 
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nondimensionals, that is, gj- ^ = 


= F , is the square of a Froude 


g & 

number . In place of © and Ap/p* 0 and F can he chosen as 
independent variables, in which case 


G 2 ftp *> ®> F ) 


(4a) 


replaces (4) . 

The simultaneous addition of two new independent variables 
resulting from the stratification means an unusual increase in the 
calculation, which in the homogeneous case is already quite extensive- 
This function G 2 is calculated explicitly numerically but not 
analytically . To arrive at a relation that is conveniently comparable 
with experiment, the disturbance wave length X Q is, ultimately, 
eliminated, by searching for that Reynolds number R = R k (critical 
Reynolds number) where precisely a single undamped disturbance still 
exists, while all others are damped. 

g 3 ( R k> ©> F ) =° 


or 


Rk = Rk(e, F) (6) 

is then also determined and from it a critical Richardson number 0 k of 
the profile, that Is, the maximum value of 0 at which an undamped 
disturbance wave length still exists at all. This depends therefore 
also on the Froude number F: ©^ = 0 k (F) - According to the energy 

considerations by Prandtl and Richardson, © k = 1 is independent of 

the velocity profilej according to Taylor and Goldstein, who disregarded 
the inertia effect of the stratification with respect to the gravity 
effect, that is, assumed F = 0, ©^ = 1/4 for certain linear velocity 
profiles and uniform density distributions. 

Tollmien's value of the critical Reynolds number for plate flow in 
homogeneous fluid was R^ = (XV 5 */ v ^k = 575> hence Rfc(0, F) = 575 
for 0 = 0 and all F. It is to be expected that (0 > O)Rk(0) > 575 for 

stable stratification, and (0 < O)R k (0) < 575 for unstable stratification. 
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The stability study with stratification of densities is, first, made 
rather general for a profile with curvature different from zero and then 
the laminar flow along a flat plate is calculated as example, so as to 
tie in with Tollmien and with experimental data- One more important 
assumption is made, namely, that the stratification exists only in the 
boundary layer, while outside of it, where the velocity of the laminar 
flow is constant, the density itself is assumed to be constant. For the 
distribution of the density, an exponential law 

<k 

p 0 (y) = p 0 w e " 7y 


hence that — = -7 is a constant, 

p iy 

mathematical simplification. 


is assumed later for reasons of 


2. THE GENERAL EQUATION OF DISTURBANCE 


Suppose the undisturbed plane laminar flow has the direction of the 
horizontal x axis and is given solely as function of the height y: U = U(y) 

(y axis at right angles, upward). The density of the undisturbed flow at 
height y is unknown: p Q = p D (y) • Limited to the two-dimensional case, 

the equations of motion and of continuity read 


/du du du\ _ dp / d 2 u d^u \ 
’U ^ = “dx + + ay2 J 


/dv dv dv\ dp / d^v d^v\ \ 

p — + u— + v—) = + p + - pg > 

\dt dx d y) dy U 2 ay 2 


du dv 
— + — = 0 
dx dy 


(7a) 


J 


The presumed incompressibility is expressed by the fact that every 
particle maintains its density during the motion, hence 


Dp _ dg 
Dt ~ dt 


& 

dx 


dp 

+ ur^ + 



(7b) 
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The plane disturbance flow Is assumed as a wave motion advancing 
in x direction^- 

We put therefore 


u = U(y) + u* 


v = v ' 


A 


where 


> (8a) 


U 1 = ui(y)e 1 ( ajc "* 3t ) v 1 = v^yOe 1 ^"^ _J 


are the components of the disturbance motion; a = 2«/x is real and 
denotes the spatial angular frequency; 0 is, in general, complex; 
p = p r + iP^. P r is the angular frequency of the disturbance motion 
with respect to time; P-^ indicates the amplification or damping, 
depending upon whether positive or negative- 

The pressure p and the density p are expressed as 


with 


p = p 0 (y) + p'. 


p = p 0 (y) + p' 



p' - 


/ \ i(ctx-pt) 

= Pi(y)e 


(8b) 


■'’This implies no limitation for the general character of the 
stability study, since Squire (Proc- Boy- Soc- A- vol- 142, 1933) indi- 
cated that when a plane flow is unstable against three-dimensional 
disturbances at a certain Reynolds number, it is unstable at an even 
lower Reynolds number for two-dimensional disturbances • The two- 
dimensional disturbances are therefore more 'dangerous" than the three- 
dimensional- 
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Inserting (8a), (8b), in (7&)> (7b) gives 

— \ 

-iapp + p ( -a 2 Up + up") 

Pi' + p(-a 2 vp + Vp") 

> (9) 

i(aU - ^)pl + V 1 P 0 * =0 —> 


Po[I(aU - P)up + U’v^j = 
p 0 [i(aU " 0)vp + gp^J = - 
Vj 1 + iaup = 0 


only terms of the first order "being maintained; the dash denotes the 
differentiation with respect to y • There are four equations for the 
four unknowns u^, v-^, p-^, and p-^, which, after elimination of up, pp, 

and pp leaves one equation for vp namely 


(aU - 3)^(vp" - a 2 vp) - a(aU - 0)u"vp 1 + — ■ T-ga 2 Vp 

Po L 

+ (aU - P) 2 vp' - a(aU - 0 ) U 1 v-^ = -iV(aU - 0)^Vp"" - 2a 2 Vj" 


* aK i + ?y(V" ' a V)| 


( 10 ) 


V = p/p 0 is assumed constant 2 . 


A flow function of the disturbance motion is introduced 


t(x,y, t) = cp(y)e i ^ a}C ^ 


p 

According to it, the viscosity p itself varies with y, and in 

p (y) dp 

such a way that — ; — r = const. But the terms with — are disregarded 

p 0 (y) • 

as small of higher order. 
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so that, as a result of 


u = 


^y 


and. v 1 



^1 = q>'(y) v-L = -iacp(y) 


Putting further c = p/a, where the real part of c signifies the" 
phase velocity of the disturbance motion, and introducing nondimensional 
variables by referring all velocities to the maximum velocity of 

the laminar flow and all lengths to the boundary-layer thickness 6 of 
the laminar flow, the- fundamental differential equation for the flow 
function cp of the disturbance motion follows by (10) as 


(U - c) 2 (cp" - a 2 cp) - (U - c)u"cp ^cp 

Um 


+ k 57 [< u - - 0 > D ’'p] = - 4 (d - <=>[»"" - 2 “V' 


^ + k ^ 2(cp "' ‘ a2<p 'i] 


+ <r 


(11) 


The dash indicates the differentiation with respect to -gj U stands 
U c 

for —, c for — , a for a6, and R = — — is the Reynolds number- 
421 I'm v 


The boundary conditions are: disappearance of both disturbance 

components for y = 0 and y = 00 , when a laminar flow is assumed that 
at y = 0 is bounded by a fixed wall, while being infinite upwardj 
hence cp = cp 1 =0 for y = 0 and y = 00 . With these boundary condi- 
tions, the stability problem proves to be a characteristic value problem 


of the kind that each specified group of four values 

6 d P Q 


a, R, 


— 

p o dy Urn 2 ' 


Po dy 


gives a complex value of c, whose imaginary part decides between 
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amplification or damping of the particular disturbance in the particular 
flow, which is characterized by a Reynolds, Richardson, and Frou.de 
number. However, the study is restricted, as stated in the introduction, 
to disturbances lying between stability and instability, hence for 
which c is purely real. For this case, c< U m and a « R, provided 
that u" does not change signs in the laminar flow- The laminar flow 
contains, therefore, a layer in which the phase velocity of the disturbance 
motion is equal to the primary flow velocity. This layer is called 
critical layer (y = y cr it) and plays a prominent part in the subsequent 
study. 


3. THE FRICTIONLESS DISTURBANCE EQUATION 


Because aR is very great at the stability limit, information 
about the solutions of the general disturbance equation (ll) can be 
obtained by analyzing the frictionless disturbance differential equation 


(U - c) 2 (q>" - a 2 cp) - (U - c)u"cp - ^ ^ ^cp 

+ ” ^[(U - c)V " (U " c)U'<p] = 0 (12) 

Visualize the speed of the laminar flow near the critical layer where 
U = c expanded in a series 


U " c = U k ’ (y - y k ) + rrU k "(y - y k ) 2 
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and "broken off with, the quadratic term. By (12) 


1 U k 

(y - yk) + p T A(y - y ^ 2 

d u k 


Uv 


"(cp" - a 2 <p) - — -(y - y k ) 

Uv 


1 u k" 


9 


where 


+ Kcp + L 


(y - yk) + §■ —,(y - yri 2 

2 V 




Uv 


i + — (y - y*> 


9 


1 ^ 

(y - yk) + - T -rr(y - yk) 2 

^ U k 


4 


= o 


(13) 


K 


_b_ d-Pp r5 1 _ ^Pq / / dU\ 2 

Pody U m 2 u k ' 2 _ Po^/W/k 




L = 


-JL ^£c 

p 0 ^y 


(14a) 


To make K and L constant in the boundary layer, the simple exponential 
distribution law 


Po = p o e 


-7y 


(15) 


is assumed for the density. Then 
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The frictionless differential equation contains the two nondlmen- 
sionals K and. L as a result of the stratification of densities, K 
representing the gravity effect and L the inertia effect of the strati- 
fication on the disturbance motion- Taylor and Goldstein ignored the 
inertia effect with respect to the gravity effect almost altogether which, 
according to (13) , seems to be permissible only when L« K, or 


according to (l4b) when 


/dU\ 2 6 V; 
W/v g ~ 8 g 


« 1, that is, when the Froude 


number of the flow is small with respect to unity. But, from the 
subsequent calculations, it will be clear that at equal numerical value 
of K and L the influence of the gravity effect on the stability is 
far greater than that of the inertia effect- „For the measurements with 


which the calculation is 


to be compared, 

8 g 


even exceeds unity, hence 


L > K- Since it is impossible to give an estimate of the amount of the 
inertia effect involved, the complete calculation is carried out with 
gravity and inertia effect- 


The next step is the integration of ( 13 ) • The solution of this 
differential equation is visualized as expansion in powers of K and L 


cp = cp® + Kcp® + Lcp® + . . . (16) 

broken off with the linear terms in K and L- It implies, proceeding 
from the homogeneous fluid, the study of a slightly stratified flow 

(K = L = 0 is the homogeneous solution) j is the already known 

solution of the differential equation of disturbance of the homogeneous 

fluid and cp® and cp® are additional solutions due to stratifica- 
tion. This expansion, and particularly the stopping with the linear 
terms, is permissible only for values of K and L which are small 
with respect to unity. 

As the laminar flow, whose stability is to be studied, is approxi- 
mated by a constant, a linear, and a quadratic function, we integrate in 
the following equation ( 13 ) for constant, linear, and quadratic velocity 

distribution. The calculations for cp® and cp® run parallel. 

For the zone of constant laminar -flow velocity, the density is 
assumed constant, as stated in the introduction. The disturbance 
differential equation then is 

cp - a5cp = 0 ' 
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with the decaying solution 

q> = e‘ a y (17) 

for great y. 

For linear velocity distribution with the new variable, it is 


and with 



ay k = a x 



(18) 


by (13) 


1 W 


a-l 2 cp] 


+ Kqp + Ly k |y x cp - y-j/ 



(19) 


For the homogeneous fluid, there follows from 

d^cp(o) 


dyn 


- a3_2q)(°) = q 


the solutions 


CPl 


( 0 ) _ 


sinh a-jy^ 


a 


<P 2 


(0) = 


cosh a. 


1^1 


( 20 ) 


By (l6) and (19) the linear inhomogeneous differential equation for 


the additional solution for cp 


(K) 


is 


jf&r ' “i 2 ” 00 ) “ ' (p<0) 


(21a) 
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and for the additional solution 9 


(L) 


yi 


/ d g cp (L) 
v 4 yi 2 


- ai g cp (L) ] = W9C0) _ yi 4 9 


( 0 )> 




( 21 b) 


the homogeneous part of both equations being the same as the differential 

equation for cp(°). Hence, only a particular integral of (21a) and (21b) 

has to be determined to produce the general integrals cp (*0 and cp(k). 
These integrals can, however, be determined immediately by integration. 

The result is a fundamental system for the solutions cp( K ) and cp(b) 
by 21 a, 21 b) in the form 


Tv 


(K) 


= + 9v 


( 0 ) Ti(°) 1 9 2 ( 0 )qx v / 0 ) 92(0) r Jl Tl^Tv^ 

Jjo 


yk 






yk 




<iyj 


9 


(L) 


+ cpv^ 





d9 v (°A 

— r 1 



+ 9 2 





yi 


d^O)' 

<iyi 


)dy 1 ( v = 1,2) 



@> • 


The choice of the lower integration limit is immaterial j the 
present choice is the point y = y Qj where the linear velocity profile 
is Joined to the parabolic. 


For parabolic velocity distribution, which is written in the form 


U = (a - y)2 


(23) 
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(a = given constant), there results from the new variables 

y - 


y 2 




and 


a g = (a - y k )a 

from equation (13). the disturbance differential equation 
^ ' ?) ^ + 72 ( x ‘ + ^ 


+ L(a - yk) 


^ " f) (1 " 72)9 ' y 2 2 ( X - 


dcp 

ayv 


= 0 


(24) 


The solution cp(^) D f this equation for the homogeneous fluid 
(K = L = 0) has already been computed by Tollmien and others. A 
fundamental system is given by 


9l^ o 

— 37 = d 1^2 + d 2^2 + ^2 


3 + . . 


91 


( 0 ) 


9< 


,(°) = e Q + ei y 2 + e 2 y 2 2 + . . . - - log y 2 , 


CXg- 


Og^ 


where dq = 1; ^ ; *3 = -g-j = ^ 


n 2n(n - l) 


n(n - 3)^-! + 2Q2 2 d^ - Og 2 d^ 


°Cg2 q Og2 

e o = 1; e l = 0; e 2 = — - 1; e 3 = g + -jgj . . 


n = 2,3, 


(25) 


3 n = gn ( n 1 - "l) ^ n ' 3) e n -l + 2a 2 2e n-2 ' °2 2e n-3 " (2n " 3)^-1 +2(2n - l)dj 


n = 2,3 • • 


NACA TM 1262 


21 


By (16) and (24), the differential equations for the additional solu- 
tions cp (K) and are 


<■ - f - ^ ,o ) * 4 - a* 


>< K ) = -cp(O) 


( 1_ f)Pt ~( a _ yk)l( 1- • ^2 


/ y 2 \dcp(°)| 
\ 2/ dy 2 J 


The homogeneous part of this inhomogeneous differential equation 
for cp (K) and <p(^) again agrees with the differential equation 
for cp(^). hence a fundamental system for the general solutions ep(®“) 
and is obtained again by quadrature 


,(K) = rrv. (0) . 


■^4 4Fi 


(0) I 2 Cp q^°)cp v ( 0 ) 

a FTW 


dy 2 (V = 1,2) 


(1 - y2)q>v^ 0 ^ (0) 

- cp,,' 7 

4'?) 


»«>-«(« „(« ', - pk a. 


+ ®J°> 


(0) CW^_ (0) , 

' r B /i - ? 
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Herewith the frictionless solutions of the stratified fluid for 
linear and parabolic velocity distributions are obtained. Of these 

integrals only cp-^^ in the singular point y = y v. of the differential 

equation is regular, while and <P k ^) , and cp 0 (k) 

exhibit singularities, where either they themselves or their derivatives 

contain terms with log(y - y k ) or (y - y k ) These singularities are 

obviously attributable to the disregarded friction, and disappear as soon 
as they are properly allowed for. The integration along the real y axis 
must be replaced near the singular point y = y k by the integration in 

the complex along a semicircle around point y = y k . First, the sense of 

rotation about point y = y k must be decided upon, that is, the choice 

of the branch of the log at transition from positive to negative y - y v . 

( 0 ) * 

Tollmien indicated this transitional substitution for cp 2 which is 
obtainable when the entire equation (11) is discussed in close vicinity 
of the critical point y = y k . 

BEHAVIOR OF THE SOLUTIONS IN THE NEIGHBORHOOD OF 
THE CRITICAL LAYER 


To this end, a small zone around point y = Yfc (transitional zone) 
is assumed, in which the substitution of U - c by U k '(y - y k ) and 

U' and u" by U k ' and U^" is sufficiently accurate. In addition 

y - y k = (oRUk'r 1 / 3 n = eq ( 27 ) 

where, on account of the small value of e, even in the transitional 
zone (that is, small values of y - y k ), q can assume great values. 

By (11) and (27), the differential equation for cp(q) is 

iqcp"" + TjCp " ( tj - 2iae^) - qcp/e-^- + e^a^T] - ie^aM 

lv 1 

+ Kcp + —— €T|(cp - T|Cp ' ) =0 

V 2 


(28) 
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ft 



1 5 d Po 

where the terms with — — 5 — in (11) are disregarded as small of the 

n p Q dy 

second order. Tollmieii (reference 13) demonstrated by means of the 
corresponding differential equation for homogeneous fluid (K = L = 0) 

that the transitional substitution for cp 2 (®) reads as follows 


7v >0 : 9p 


( 0 ) = 


+ e 


iy + e 2y c 


y - yfe <0:cp 2 (°) = e Q + e x y + 


e 2 y 2 + 


Th 

Ub 


U- 


- 5 qp 1 C°)log(y - y k ) 


+ T *r» cp i^ 0 ^ log |y ■ - iir ) 

U-i_ 


(29) 


Thus, of the infinitely many values of log(y - y k ) + 2krti 
(k = positive or negative whole number), the one with k = -1 proves 

physically real. A calculation similar to that made by Tollmien for cp 2 ^^ 
gives the same result for cp^v^) , cpg( K ) , and cp 2 ( L ) , namely, that the 
termB with log (y - y k ) for negative y = y k must be replaced by 
logj y - y k j- in. .A further approximate calculation for the homogeneous 
fluid, which in (28) includes only the highest friction term (icp"" + qcp" = 0), 
gives the two other solutions cp^ 0 ) needed for representing 

the general integral of the differential equation of disturbance of the 
fourth order (ll). It is 



, 1 % 5 ,1W2) 




dq 


(30) 


hC 1 ) »( 2 ) is the Hankel function of the first and second kind, respec- 
tively. 

Additional solutions for cp 3 and 94> namely. 


93,4 = 93,4(0) + Kcp 3 ^( K ) + i/p 3jk ( L ) + . . . 

should also be computed, but the subsequent paragraphs indicate that 
the calculation of 93,4^ and 93,4^^ is superfluous because 




I 
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H,2 (K) 

, <P 3 ,l/ L) _ 91,2 

and — « — i — 


^ 2 (0) 

”3,4 (0) ” 1,2 


(L) 

( 0 ) 


(31) 


The effect of the stratification of densities on the frictionless 
solutions <P]_ g therefore, substantially greater than on the friction 

solutions cp^ 4* so ‘the latter can be disregarded in the approximation; 

this simplifies the calculation considerably. The result is analagous to 
the conditions in a flow with stratification of centrifugal forces. In 
the earlier Gottingen report (reference 20), it also had been shown that 
in stable stratification due to centrifugal forces the effect of this 
stratification on the frictionless solutions 9-^,2 very much 

greater than on the friction solutions qp^ 4 . 

The following proof is given: 


5 • THE VARYING EFFECT OF THE STRATIFICATION OF DENSITIES 
ON THE TWO PAIRS OF SOLUTIONS 


Desired is a representation of the four integrals cp^, cpp> T3, 
and 94 for the wall proximity (very small y) on the specific assumption 
that the phase velocity c of the disturbance motion is very much higher 
than the velocity near the wall . In other words, the critical point is 
to lie outside the layer where the solutions are analyzed . For the wall 
proximity U - c can then be replaced by -c and U by U w . 


1 dp o 

On account of 5 — = 7, according to equation (15), the general 

P Q a Y 

equation of disturbance for a layer near the wall is by equation (11) 


cp"" + (iaRc - 2a^)cp" + + iccR^Uw" + ~ 9 = 0 

The solutions of this equation with constant coefficients are of the form 




V 


= e k vy 


v = 1, 2, 3, V 
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k v are the four roots of the equation of the fourth degree 

k^ + (iaRc - 2a2)k 2 + ja^ + laR ^U w " + ~ ct^cjj = 0 


Both pairs of solutions are given "by 


k 2 = -i(ic®c - 2a 2 ) ± ^j-(iaRc 2a 2 ) 2 - - iaR^U w " + p - 


2 

a c 


hence 


*1,2 -“ 2 


+ ?) 


or when XL, = 0 
w 


k l,2 2 " a2 


1 - E 


(5 l 

\ac 


i\2 


, since K = 


7K 

rr .2 
u w 


or 


k l,2 = 4 a 


1 - 


K/U. 


t \ 2 




ac 


Likewise, when a 2 is disregarded relative to aRc 


2 /V 

k 3,4 = - lot:Rc + \~T 


or 


!,4 = + \J- 


IaRc 


. Kia/VV 

2Rc\ac / 


From 


cpi = + Kcp]_( k ) = 


a H - — 

2 \ac 


Kf U w 


• \ 2 
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follows 


*1 (K) oA’Y 


CPl 


( 0 ) 


2 \ac , 


and likewise from 


cp^ = «p (°) + Kip- 


(K) _ 


\JictRc 


= e 


Kf a\/ U w ,N 


2 \Rc/\clc 


J 


c P 3 (K) _ + /ia 

q, (0) " + 2\oc / \/ Rc 


Hence 


(k) cp. Con! 


93 ''-’ Tq 


X 


q>i^ cp 3 (°) 



(32) 


With 6 as reference length, the flow along the flat plate 1 b, at the 
stability limit, (reference 22) approximately 


a = 0.8V; R = 1700; c = O.V-2 


hence 


<p M cpJO) 

X 


cp 1 ( K ) Cp^®) 


— ~ — « 1 
850 29 


The same holds true for and > and a similar derivation 

can he given for and cp^^) , with which equation ( 31 ) is therefore 

proved . 
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m 

// 


This indicates that the effect of the stratification of densities 

±cl y 

on the slowly variable solutions cp^ and qp 2 , which are as e , I s 
much greater than on the rapidly varying solutions which are as e ± V'aRcy . 


6. FORMULATION OF THE CHARACTERISTIC VALUE FROBLFM 


The next problem deals with the "boundary conditions and also the 
formulation of the characteristic value problem. A basic velocity 
profile is assumed which Increases from value zero to a maximum value U^ 
which it then maintains constant. By (17), the solution of the 
differential equation of disturbance for the zone of constant velocity, 
where constant density is assumed also, is: cp = e _a, y. With y = a 

denoting the point of contact with the zone of constant U, the boundary 
condition at this point reads 

V + ac P a = 0 (33) 

The general integral reads 

<P = Cq'Pq + C 2 (p 2 + C 3 <p 3 + C^q > 4 

The boundary conditions (cp = cp 1 =0 for y = 0 and y = 00) can be 
considerably simplified on the basis of the particular properties of 
the four integrals cp 1 , cp~, cp^, and cp^. Thus, the integral cp^, 

^ — 1 

which for very great y increases Infinitely (like e ’ 07 ) need 

not be considered in the general integral, so that C4 = 0. Furthermore, 

the integral CP3, which Is as e needs to be considered only at 

the wall, while In the connecting point on the zone with constant velocity 
(y = a) it can be regarded as zero with good approximation. 

For y = a, the boundary conditions read, by (33) therefore, simply 


Cl'P’la + c 2?2a + a < C l < Pla. + C 2^ = 0 


or, when putting 


?va + ^a =J Va 


v = 1,2 


(33a) 
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C-,#. +• C 

1 la 


2 $ 2a ~ 0 


(34a) 


The "boundary condition at the wall (y = 0) reads 


C l^lw + c 2 c P2w + C 3 c f ) 3w - 0 

(34b) 

C l cp 'lw + C 2 (p '2w + C 3 Cp '3w = 0 

(34c) 


hence the three homogeneous equations (34a), (34b), and (34c) as boundary 
conditions. If one solution qp other than zero exists, the dete rminan t 


•la 

$ 2a 

0 



'Plv 


<P 3 w 

= 0 

(35) 

'P'lw 

*'2w 

cp ' 

Y 3w 




must disappear, wherewith the stability investigation as a characteristic 
value problem of the differential equation of disturbance is proved. 

The complex equation (35) depends, besides the constants of the profile 
used as basis, upon the parameters a, c, B, K, and L - all of which, 
especially c, are purely real - since only the transitional points from 
stability to instability are involved. The functions <p Vw and 3> Va are 
usually complex (compare equation (29) and (30))- Visualizing the 
parameter c eliminated from the two real equations with which the 
complex equation (35) is equivalent leaves one equation in which a, B, K, 
and L are contained. This equation gives, for fixed K and L, the 
indifference curve in the oE plane corresponding to this stratification 
(K, L), and which separates the damped from the undamped disturbances. 

Thus the solution of the stability problem terminates in the discussion 
of equation (35)* 

Calculation of the determinant (35) gives 

c f , 3w _ ^2w^la ^lw^a 
^ ' *'2^1a - <P V 2a 


( 36 ) 
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This equation is then expanded in powers of the stratification 
quantities K and L corresponding to the fact that, starting from the 
homogeneous fluid, the stability study for slight stratification is to 
be carried out. Only the right-hand side of the equation has to be 
expanded, which contains the frictionless solutions exclusively, hence 
is not dependent on the Reynolds number R. On the other hand, the 
solution cp^ no ^ dependent on the laminar flow serving as basis. 

It can be summarily taken over from Tiet Jens 1 report . The left-hand 
side of equation (36) depends only on 

■Hw = ™ (37) 


For D(t} w ) = - — 
e 

Tiet Jens: 



the following table I is obtained according to 


Table 1: D(q w ) 


_T lw 

D(^lw) 

0 

0.T02 - 0.425 i 

0.5 

O.785 - 0.1+11 i 

1.0 

0.920 - O.389 i 

1.5 

1.043 - 0.297 i 

2.0 

1.206 - 0.147 i 

2.5 

1. 357 + 0.108 i 


as Function of 


-7W 

D(tJv) 

3.0 

1.400 + 0.515 i 

3.5 

I.180 + I.130 I 

4.0 

0.460 + 1.250 i 

4.5 

-0 . 0405 + O.8080 1 

5.0 

0.0057 + 0.3645 i 

5.5 

0.1913 + 0.2393 i 


After putting 


1 TW^la " <P 1 w <I, 2£ 


7 k q>' ■ #. 

2w la 


ro ' <I> 

^ lw 2a 


= E(a, c, K, L) 


equation (36) gives 


PQjy) 

^W 


E(a, c, K, L) 


(38) 


(39) 
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This complex equation forms the starting point for the subsequent calcu- 
lations, the aim of which is the explicit representation of the relation 
"between the five quantities a, c, R, K, and L given by this equation. 


7 • SOLUTION OF THE CHARACTERISTIC VALUE PROBLEM 


For the numerical treatment of this equation, it is essential that, 
of the five parameters on the left-hand side, the two stratification 

quantities K and L do not occur nor the Reynolds number R on the 

right-hand side . 

The solution of this equation, that is, the calculation of the 
correlated characteristic values a, c, R, K, and L, is first effected 

analytically so far as the expansion in powers of K and L is 

involved. The subsequent treatment of the ensuing equation must be 
made by numerical graphical method, like in the earlier stability 
studies (references 21, 22) . Since the whole calculation in its 
analytical, as in its numerical-graphical part, is far too extensive to 
be reproduced here, a brief outline of the line of reasoning must be 
sufficient : 


The expansion of equation ( 39 ) in powers of K and L is formally 


^(Hw) \ 

E^Cct^c) 

Hw 


+ 



+ 



+ 


(39a) 


E^°) denoting the value of equation ( 38 ), when inserting for tp^, cp^ 
the solutions for the homogeneous fluid, that is 


E (°> = -i- 

yk 


(0L (0) _ (OK (0) 


( 6) 7 

*2w *l £ 


To T 


'Plw' 


"2a 


CPl- 


(0)K (0) 


w 


"2 a 


(40) 


cSE dE 

while — and — are obtained when introducing for all terms 
cK SL 

their expansion in powers of K and L, that is 


in (38) 
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<Pw = Cp vw ( ° ) + Kc Pvw (K) + Lt Pvw (L) + 

®Va = ®va^ 0) + K 3 > Va (K) + L 4 va (L) + 



= 1 , 2 ) 


(41) 


q> v (^), and hence ® and ® V a^^ are ^ movn according to 

section 3* 


The complex equation (39) Is equivalent to two real equations. For 
the left-hand side of equation (39), the decomposition into real and 
imaginary part is immediately given according to table 1. The right-hand 
side contains an imaginary component as a result of the transitional 
substitution of the integral q >2 (equation (29))* Effecting the indicated 
expansion in powers off K and L explicity results, after considerable 
paper work, in the expressions for the real and imaginary part of 


E 


5, «* 2 

3k 3l 


numerical values of 


as function of the two variables a and c . The 
r n \ 3e 3E 

E^ u ' , — , — , for the velocity profile with which 


3k 3l 


the calculation was made, are indicated in table 2. The further solution 
of equation (39) is made by graphical method. With fixed c, K, L the 
left-hand side of equation (39) is plotted with q w as parameter, and 
the right-hand side with a as parameter in a polar diagram (fig. 2) . 

In general, it results in two intersection points of the — curve with 

^w 

the E(a, c, K, L) curve for which the parameters a^, and 

are found. The corresponding Reynolds number follows then from the equation 


n v = -yk(aRUk') 1 ^ 3 


The great advantage of the expansion in powers of K and L is 
manifested in the numerical calculation by the fact that the relation 
of the desired transition points from stability to instability of the 
stratification condition is verified in very simple manner. Computing 

/ 0 x 

the quantities E' u ', — , — for fixed c, as function of a, immediately 

cK oL 

yields the curves E = E(a,c,K,L) for the homogeneous flow and for all 

rW Vtji 

states of stratification. From the numerical values of ^ and ^ in 
table 2, it is apparent that the influence of the gravity effect of the 
stratification is very’ much greater than the influence of the inertia 


effect 


/se aErt 
\5k >:> 3L /• 


In the present calculations, small negative values 


of K and L were also included, that is, slightly unstable stratifica- 
tions, such as the flow of an evenly tempered air current above a heated 
plate, which is a little warmer than the air. 



uo 

ro 


TABLE 2.- R(E), l(E), r(|§), l[g|), r(||), i(j|)aS FUNCTION OF a 
FOR c = 0.20; 0.25’; 0.30; 0.35; 0.375; O.to; 0.425; 0.45. 


a 


R(U) 

1(E) 

~W 

t / 5 e\ 

J8e\ 

w) 

<£) 


R(E) 

1(E) 

Jl l) 
R ml 

t (Ze) 

i \ 5 k ) 

Jte\ 

v§rj 

t/Se) 

HSl] 





c = 

).20 




c = 0.25 

0.20 


0.615 

0.009 

10.00 

2.03 

0.820 

0.021 


0.546 

0.026 

8.64 

2.11 

0.730 

0.010 

.30 


.396 

.024 

10.88 

3.22 

l.4i 

.050 


• 350 

.050 

9.44 

2.80 

1.17 

.000 

.to 


.132 

.0V7 

11.78 

k. to 

2.10 

.095 


.141 

.090 

10.28 

3-59 

1.49 

-.005 

.50 


-.205 

.085 

12.8^ 

5.94 

3 .o 4 

.159 


-.165 

.155 

11.46 

4.37 

2.18 

-.oto 

.60 

.70 


-.610 

.150 

14.06 

7.79 

4.10 

.263 


-.513 

.273 

12.56 

5.08 

2.89 

-.155 



c = 0.30 


c = 0.35 

0.30 


0.642 

0.026 

7.10 

1.32 

0.490 

-0.034 


0.730 

0.027 

5.92 

0.850 



.to 


.toi 

.052 

7.74 

1.77 

.690 

-.072 


.615 

.052 

6.45 

1.10 

o.46o 

-0.110 

•50 


.333 

.091 

8.38 

2.15 

.800 

-.135 


.480 

.090 

6.96 

1.28 

• 570 

-.192 

.60 


.128 

.157 

9.28 

2.55 

1.100 

-.255 


.335 

.151 

7.56 

1.47 

.7to 

-.335 

.70 


-.111 

.260 

10.26 

2.87 

i.46o 

-.465 


.165 

.245 

8.30 

1.56 

.910 

-.565 

.80 


-.too 

.kk2 

11.80 

3.84 

1.820 

-1.09 


-.029 

.402 

9.08 

1.52 

1.080 

-.985 



c = 0.375 


c = o.to 

o.to 


0.652 

0.050 

5.89 

0.820 

0.369 

-0.118 


0.690 

0.050 

5.48 

0.586 

0.324 

-0.123 

.50 


.5to 

.088 

6.33 

• 925 

.500 

-.208 


.590 

.087 

5.76 

.650 

.430 

-.219 

.60 


,to7 

.ito 

6.85 

• 939 

.632 

-.350 


.484 

.139 

6.14 

.645 

.520 

-.355 

.70 


.285 

.232 

7.32 

.966 

.742 

-.573 


.371 

.222 

6.58 

.506 

.590 

-.570 

.80 


.133 

.369 

7.91 

.803 

.812 

-.923 


.250 

• 339 

6.91 

.346 

.595 

-.860 

.90 


-.024 

.574 

8. to 

.340 

.870 

-1.402 


.129 

• 525 

7.26 

.056 

.452 

-1.26 



c = 0.to5 


c = 0.45 

o.to 


0.720 

0.0 to 

4.96 

0.380 

0.270 

-0.125 


0.754 

0.046 

4.58 

0.170 

0.237 

-0.127 

.50 


.637 

.085 

5.2k 

.390 

.350 

-.223 


.681 

.082 

4.84 

.ito 

.288 

-.220 

.60 




5.51 

.330 

.423 

-.355 

1 

.605 

.126 

5.00 

.070 

.336 

-.342 

.70 



.209 

5-94 

.190 

.44o 

-.555 

I 

■ 530 

.195 

5-36 

-.050 

.324 

-.513 

.80 


.360 

.315 

6.12 

.000 

.395 

-.795 

I 

« 456 

.298 

5. to 

-.030 

.222 

-.725 

.90 


.266 

•to5 

6.28 

-.220 

.210 

-1.10 


.391 

.430 

5.44 

-.490 

.014 

-.927 
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8. APPLICATION TO PLATE PLOW 


For the stability study with stratification. of densities, the 
laminar flow along a flat plate has "been chosen as example for the 
numerical calculation. Tollmien made the corresponding investigation 
for the case of homogeneous fluid (reference 13 ) • Blasius originally 
computed the laminar velocity distribution along the flat plate 
(reference 23) according to formulas by Prandtl (reference 2h) . In the 
present stability study, the modification in the laminar flow due to 
the stratification of densities can be disregarded because a simple 
evaluation^ indicates that it is small, when the boundary layer contains 
only a slight stratification of densities as it was assumed. 

- 

The force in x direction produced by the stratification per unit 
volume is: 


(pg) 


1 dp R dS _ dS 

0 aF s 5 ' p ®' 6 E 


and the friction per unit area n— 

dy 


or per unit volume 


M. dU 
6 dy 


PU m 




Accordingly, the effect of the stratification of densities on the 
laminar profile can be disregarded when 


P87° — « 

dx 6^ 


In that case 



hence it muBt be 


« ^ Um 


or 


gy_ _ 


«1 


5 


dB 

dx 


v 

Ik 


which is precisely the assumption made above. 
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For the present purposes, an approximate representation of the 
Blasius profile "by a straight line and a parabola is used (fig. 1). 
The profile is to start at the wall (y = 0) with a straight piece, 
which is tangentially joined to a parabolic piece and it, in turn, 
tangentially onto the constant velocity U = U m . 


The result is the approximate representation 


0 < | < 0.175 ~ = 1.68 | 

6 U m 5 


0.175 < | < 1.015 ^ = 1 - (l. 


015 


(42) 


z = 1.015 JL = 1 


u, 


in 


earlier indicated for U. 


Hence a = 1.015 by comparison with equation (23). But for the 
transitional substitution (equation (29)), U"/U’ must be taken according 
to the exact Blasius series, not U' = constant and U" =0 in the 
linear range and U" = constant in the parabolic range. The first terms 
give 


Ul 

u m 


1.68 1 


E 



IT 

U, 


m 


-18. 



(^3) 


This supplies all the data necessary for the numerical calculation. 
Figure 2 shows the previously described polar diagram for a specific 
c-value (c = 0.10), where the several E curves correspond to the 
written-in values of the stratification quantity K, and L is always 

equal to zero. For this c-value - say, for K >— r- - there are no inter- 

% T) 

section points of the E curves with the _ curve, hence, no 

'Hw 

instability. For other c values, this maximum K value, at which 
instability is still just possible, is different. The Highest possible 
K value at which instability is still just possible (critical value 
of K) is K = l/24, when assuming L = 0. For other values of 
L (L = -l/lO, +l/l0, +l/5), other critical values of K are obtained. 
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c 

1 : 

(aS*)l 

0.20 

1.15 X 10 ^ 

0.038 

0.30 

1.84 x 10 3 

0.078 

0.35 

1.01 x 10 3 

0.101 

0.375 

7.25 x 10 2 

0.121 

o.4o 

5.86 X 10 2 

0.139 

0.45 

4.o4 x 10 2 

0.193 



K = ~S5 

0.20 

0.30 

0.35 

0.375 

0.40 

0.45 

9. 27 x 10 3 
1. 68 x 10 3 
9.20 X 10 2 
7. 15 X 10 2 
5.70 X 10 2 
4. 51 x 10 2 

0.050 

0.093 

0.121 

0.140 

0-159 

0.233 



8.75 x lo3 0.187 
2.99 x 10 3 0.232 
1.88 x 10 3 0.251 
1.30 x 10 3 0.266 
6.08 x 10 3 0.282 


8.22 

X 

10 3 

1.60 

X 

10 3 

8.96 

X 

10 2 

7.05 

X 

10 2 

5.76 

X 

10 2 

5.18 

X 

10 2 

7.20 

X 

10 3 

3.01 

X 

10 3 

1.53 

X 

10 3 

1.15 

X 

10 3 

8.93 

X 

10 2 

7. 36 

X 

10 2 

6.33 

X 

10 2 

5.76 

X 

10 2 



6.75 

X 

io3 

2.61 

X 

10 3 

1.70 

X 

10 3 

1.20 

X 

10 3 

5.23 

X 

10 2 


3.76x10^ 0.149 
1.20 x 10^ 0.188 
4.64 x 10 3 0.223 
3.29 x 103 0.238 
2.07 x 10 3 O.251 
1.42 x 10 3 0.264 
1.02 X 10 3 0.274 

7.13 x 10 2 0.273 




2.24 x 105 
1.71 X 10 1 *' 
3. 14 X lo3 
1. 47 X 10 3 




0.2251 7-18 X 10 3 | 0.190 

























TABLE 3 — Continued 


u> 

ON 


K = - 


SO 


L = - 


HT 


K = 0 


c 

(r) 

(aS*)-> 


(a 6*) 2 

c 

w. 

(a5*) 1 

( w /)« 

(aS*) 2 

0.30 

0.35 

o.4o 

0.425 

0.45 

1.94 x 10 3 

9.94 X 10 2 
6.49 X 10 2 
5.20 X 10 2 
4.58 X 10 2 

0.072 

0.102 

0.130 

0.155 

0.184 

3.16 X 10 3 

1.38 X 10 3 
9.41 x 10 2 

6. 38 x 10 2 

0.218 

0.248 

0.254 

0.252 

0.20 

0.30 

0.375 

0.40 



8. 39 x 10 3 
1. 67 X 10 3 
7. 94 X 10 2 
6. 71 x 10 2 

0.065 

0.113 

0.166 

0.187 

5.21 X 10 3 
1.51 x 10 3 
1.07 X 10 3 

0.207 

0.250 

0.257 

K = w 


0.20 

0.25 

0.30 

0.35 

o.4o 

0.425 

0.45 

1.17 X 10^ 

4.11 X 10 3 

1.81 X 10 3 

1.03 x 103 

6.22 X 10 3 

5.29 x 103 
4.92 X 10 3 

0.040 

0.062 

0.086 

0.114 

0.150 

0.173 

0.218 

- 

7.46 x 10 3 
2.76 x 10 3 

1.28 X 10 3 

8.63 x 10 2 
5.44 X 10 2 

0.185 

0.225 

0.253 

0.262 

0.245 
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TABLE 3 — Continued 


K = lfe 


c 


(aB*)]. 

w. 

(aS*) 2 

c 

W, 

(aB*^ 


(aS*) 2 

\ 

0.20 

6. 16 x 10 3 

0.096 

1.98 X 10^ 

0.187 

0.20 

8. 53 x 10 3 

0.175 

l.l4 x 10^ 

0.199 

0.30 

1.45 x 10 3 

0.139 

4.09 X 10 3 

0.242 

0.30 

1.80 x 10 3 

0.210 

2.72 x 103 

0.257 

0.375 

6.47 x 10 2 

0.195 

1. 33 X 10 3 

0.286 

0.325 

1. 38 x 10 3 

0.228 

1.94 x 10 3 

0.265 

o.4o 

5.49 x 10 2 

0.214 

9. 75 x 10 2 

0.296 

0.35 

1. 19 x 103 

0.259 

1. 30 x 103 

0.270 

IT\ 

CVJ 

d 

4. 65 x 10 2 

0.251 

6. 65 x 10 2 

0.300 
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0.20 

9. 58 x 10 3 

0.188 

1.04 x 10^ 

0.195 

0.20 

7.12 X 10 3 

0. i4i 

1.49 x lO 4 

0.188 

0.25 

3.82 x 10 3 

0.203 

5.05 X 10 3 

0.225 

0.30 

1. 57 X 10 3 

0.183 

3.10 x 10 3 

0.256 

0.30 

1.91 x 103 

0.224 

2.58 X 10 3 

0.252 

0.375 

7. 86 x 10 2 

0.251 

1.07 x 10 3 

0.288 

0.325 

1. 55 x 10 3 

0.246 

1.76 X 10 3 

0.262 

0.39 

7. 66 x 10 2 

0.272 

8.45 x 10 2 

0.284 
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0.25 

4.10 x 10 3 

0.212 

4. 70 x 10 3 

0.224 






0.30 

2.06 x 10 3 

0.235 

2. 37 X 10 3 

0.250 
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c 


(aS*) 1 

(T\ 

(aS *) 2 

0.30 

0.325 

0.35 

0.375 

1. 77 x 10 3 
I.27 x 10 3 
9.26 X 10 2 
7. 52 x 10 2 

0.24-7 

0.256 

0.268 

0.290 

K = i 
37 

2.26 X 10 3 
1.78 X 10 3 
1.43 X 10 3 
1.04 x 103 

0.276 

0.294 

0.308 

0.319 

0.30 

0.325 

0.35 

1.88 x 10 3 
1. 35 x 10 3 
1.01 X 10 3 

0.258 

0.266 

0.280 

2.17 x 10 3 
1. 71 X 10 3 

0.280 

0.296 




K _ _1_ 

35 


0 

0 

0 



.30 2.02 x 10 3 
.325 1.40 X 10 3 
.35 1.04 X 10 3 


.2jk 

,27k 

.290 


1. 65 X 10 3 
1.26 X 10 3 
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So far, all lengths had been referred to the boundary- layer 
thickness 8 (half the width of the parabola inscribed in the Blasius 
profile, fig. l) . For the presentation of the mathematical results, the 


displacement thickness 8* = 



is chosen as reference 


length (8* = 0.3^18). 


The transition points computed for the various stratifications 
(characterized by K and L) are correlated in table 3, and the 
corresponding curves plotted in figures 3* and 6. The area 

circumscribed by the curves represents the area of unstable disturbances. 
The diagrams indicate the looked- for connection between the wave length 
of the disturbance X = 2ir/a, the Reynolds number R, and the stratifi- 
cation (K and L) as formulated in equation (4) in the introduction. 

Each diagram corresponds to a specific value of L(L = -l/lO, 0, +l/lO, 
and + 1 / 5 ), while the several curves of each diagram relate to different 
values of K. Each curve is therefore characterized by two parameters K 
and L. The curve K = 0, L = 0 is Tollmien's indifference curve for 
homogeneous fluid. For stable stratifications (K > 0, L > 0), the 
area of the unstable disturbances is smaller and for unstable disturb- 
ances greater than for homogeneous fluids. 


The indifference curves themselves exhibit the following unusual 
fact: For homogeneous fluid and for unstable stratifications each 

indifference curve extends to infinity at great Reynolds numbers, while 
for stable stratifications the indifference curves are closed. With 
increasing stable stratification, the unstable range of the wave length 
of disturbance enclosed by the indifference curve becomes consistently 
smaller, until finally at a K value dependent on L the indifference 
curve shrinks to a point located in the aR plane. 

The critical Reynolds number corresponding to a certain state of 
stratification is of particular interest. In figure J, the critical 

Reynolds number I m j is plotted against 0, for several values of 

V V 4 

u m 

the Froude number rr* Thus is the relation important for the com- 

8*g 

parison with experiment, as formulated in equation (6). For small 
values of 0, the variation of the Rv — curves with the Froude number 

o 

U m 

is rather small.. But the stability limit 0^ - that is, the 

maximum value of the stratification quantity at which turbulence is 
still possible - still depends to some extent on the Froude number. 



TABLE 4 


THE CRITICAL REYNOLDS NUMBER 


AS A JUNCTION OF THE RICHARDSON 


NUMBER 0 


0.0208 

o.oi 42 

0.0087 

0 

0.0090 

0.0153 

0.0233 

0.0314 

0.0392 

0.0409 

0.0417 


2 uJ 2 

| AND THE FROUDE NUMBER 



- 0.0200 

- 0.0100 

0 

0.0100 

0.0200 

0.0250 

0.0290 

3 crit = °-° 3 T 0 


- 0.0200 

- 0.0100 

0 

0.0100 

0.0200 

0.0250 

0.0290 

®crit = 0.0340 


- 0.0200 

- 0.0100 

0 

0.0100 
0.0200 
0.0250 
®crlt = 0.0297 


-0.0200 450 

-0.0100 500 

0 575 

0.0100 693 

0.0200 980 

0.0250 1305 

5 it = 0.0290 1850 
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The dependence 0}c = 0k(s§-) is reproduced on figure 8. For ^3- = 0 - 

\ b g / 5*g 

hence, for the case that the inertia effect of the stratification is 
disregarded relative to the gravity effect - it results in 0k = 0 .0409 • 
Therefore, the inclusion of the friction for nullifying the turbulence 
affords a considerably lower value of the stratification quantity than 
without friction, according to Taylor and Goldstein, who obtained 
^ * 0 . 250 . 


By equation (l4a) the stratification quantity K is so defined 
that the velocity gradient must be taken at the critical point which 
differs according to the position of the critical layer. But, in the 
present case, it was preferred to form the Richardson number with the 
velocity gradient at the wall, namely 


0 = 



m 


as obtained from K 



little from the K 


by the reduction 



is always close to unity, 
value s . 


the 


0 values differ 


9. COMPARISON OF THEORETICAL WITH -EXPERIMENTAL DATA 


Velocity and temperature distributions in a flow with stratification 
of densities, which are very appropriate for the comparison with the 
present theory, were made some years ago by Reichardt, in the warm-cold 
air tunnel of the Kaiser Wilhelm- Institute for flow research at Gottingen. 
The tunnel is l6 m long, 1 m wide, and 25 cm high. The velocity profiles 
were measured at the end of the tunnel with a hot-wire survey apparatus 
which, to compensate for the temperature effect, consisted of two hot 
wires of different diameter. One of these wires could at the same time 
be used for temperature recording. The study so far included only the 
stable case, where the upper plate of the tunnel was heated with Bteam 
(100°) and the bottom plate cooled by tap water (10° C), hence at a 

Not published so far in detail; see the test reports (reference 26). 
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o # 

temperature gradient of 90 C at 25 cm. Very low airspeeds* of around 
1 m/sec were indicated, since at higher speeds no perceptible effect 
of the temperature stratification of the flow is produced. Figure 9 
represents a velocity and temperature distribution measurement. In 
the central part of the tunnel where the speed is practically constant, 
the temperature is also nearly constant. 

For the comparison with the present theory, the measured profiles 
were evaluated to the effect that for each profile the maximum 
velocity U m and the boundary- layer thickness 8 (and hence the 
displacement thickness 8*) of the boundary layer are obtained. With 
these, the theoretical velocity distribution of the laminar flow was 
then ascertained and their velocity gradient at the wall 



computed. In figures 10 to 14 are shown for various maximum veloci- 
ties Uju, the measured velocity and temperature distributions for 
the boundary layers on the warm and cold plates together with the 
distributions upon which our stability calculations were based. The 
relation between temperature and density is, on account of the constant. 


^ or 

Pw T p dy T dy 

T, the absolute temperature; subscript w, the value at the wall. The 
density distribution assumed for the present theory (equation (15)) 
according to the law 


Po 


P 0 e 

w 


-77 


corresponds therefore to a temperature distribution 


T = (45) 

The constant y was determined from the measurements in such a 
way that the recorded temperature distribution agreed as closely as 
possible with the distribution given by equation (45), which served as 
basis of the present calculations (cf. figs. 10, 11, 12, 13, and l4) . 
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With it the Richardson number 


o = ZS 

w . < 2 . 2 
/dU\ * /dU\ 

W/w W/w 


g EE 

T dy 


the Froude number 


and the Reynolds number R - 


U m 6* 


for each 


5 g V 

profile can be computed. The value corresponding to the mean tempera- 
ture in the boundary layer was chosen as the value for the kinematic 
viscosity v. The numerical data are reproduced in table 5 • 


The decision of whether the recorded velocity profile was laminar 
or turbulent was made from oscillographic records of the voltage 
fluctuations in the hot wire, but it could also be observed direct from 
the shape of the velocity profile. Despite the finite tunnel width and 
the density differences, the measured velocity for the laminar profiles 
is in good agreement with the Blasius plate flow, while the turbulent 
profiles have almost the form of the 1/7 power law. 


For comparison, the test points were plotted in an R,0 diagram, 

© being expressed at logarithmic scale (fig. 15) • To each profile in 
figures 10 to 1^, there corresponds one point in this diagram, the 
laminar and turbulent profile being characterized by different symbols. 
The solid curve is the theoretical stability l imi t which gives the 
critical Reynolds number as function of the stratification quantity 0 
for the Froude number F = 0. The curves for the other Froude numbers 
are omitted since they do not differ appreciably from the curve 
for F = 0, according to figure 7. The solid curve gives the boundary 
between the theoretically stable (laminar) and the theoretically unstable, 
(turbulent) attitudes. The Taylor-Goldstein value ® cr -j^ = 0.250 and 

the value R cr it = 575 for plate flow in homogeneous fluid are also 
included. The comparison of the presented theory with the test data at 


which the Froude number 


U„ 


2 
m 
Sg 


ranges between 0.86 and 3*82 (compare 


table 5) indicates that in the theoretically stable zone only laminar 
attitudes, and in the theoretically unstable zone, with one single 
exception, only turbulent attitudes were observed.' 5 Incidentally, it 
should be noted that it is not contrary to theory when laminar states 
are still observed closely above the stability limit R cr ^^_ = 


5 

From the oscillographic records on the rejected profile, 
tinction could be made between laminar and turbulent. 


no clear dis- 


■fr 

■£- 


TABLE 5-- COMPARISON WITH TEST DATA 


Figure No. 


8cm 

8*cm 

v cm 2sec 1 

(T w “ T m)° 

J “1 

7 cm 


e 

■ 

theoretical 

experimental 

lU i 

W. 

PI. 

127.5 

4.34 

1.48 

0.202 • 

64 

0.0431 

934 

0.0173 

3.82 

stable 

laminar 

k. 

*. 

PI. 

134 

7.20 

2.45 

0.151 

-25 

0.0117 

2175 

0.0117 

2.54 

unstable 

turbulent 

13 \ 

f 

w. 

PI. 

105 

4.87 

1.66 

0.204 

60 

0.0359 

855 

0.0269 

2.31 

stable 

laminar 

L k. 

PI. 

111 

6. 13 

2.09 

0.151 

-25 

0.0138 

1536 

0.0144 

2.05 

unstable 

turbulent 

12 • 

w. 

PI. 

102.5 

5.89 

1.94 

0.204 

60 

0.0307 

975 

0.0324 

1.88 

stable 

laminar 

[k. 

PI. 

104 

6.34 

2.16 

0.155 

-32 

0 

b 

H 

® 

1450 

0.0218 

1.74 

unstable 

turbulent 

H 

V. 

PI. 

91 

6.38 

2.18 

0.204 

& 

0.0274 

973 

0.0467 

1.32 

stable 

laminar 

l k - 

PI. 

98.5 

6. 04 

2.06 

0.154 

-30 

0.0165 

1317 

0.0217 

1.64 

unstable 

partly laminar, partly turbulent 

10 \ 

V. 

PI. 

81 

7.78 

2.65 

0.204 

60 

0.0225 

1052 

0.0721 

0.86 

stable 

laminar 

[k. 

PI. 

87.5 

5.75 

1.96 

0.154 

-30 

0 

b 

H 

nos 

0.0261 

1.36 

stable 

laminar 


v. PI. = hot plate k. PI. = cold plate 
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for this curve indicates the attitudes where the amplification of the 
disturbance motion is exactly equal to zero. For an actual appearance 
of the turbulence, the small disturbances must, however, be considerably 
amplified; and that is the case only when the Reynolds number is somewhat 
greater (compare reference 22). The agreement is therefore complete. 


The value 


© 


crit 


0.250 obtained by Taylor and Goldstein is far 


from the © values denoting the boundary between laminar and turbulent 
flow according to Richardson's measurements. It is true that Taylor's 
and Goldstein* s stability investigations had been based upon different 
velocity profiles than realized here in experiment. The presented 
theory, on the other hand, whose assumptions are rather closely 
satisfied by the measurements, indicates a very satisfactory agreement 
with the experiment, according to figure 15 • 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 


46 


NACA TM 1262 


REFERENCES 


1. Richardson, L. F. : Turbulence and, Vertical Temperature Difference 

near Trees. Phil. Mag. 49, p. 8l, 1925* 

2. Taylor, G. I. : Turbulence. Joum. of the Meteoro logical Society, 

vol. 53, p. 201, 1927. 

3. Bjerknes, V., Bjerknes, J. , Solberg, H. , and Bergeron, T. : 

Physikalische Hydro dynaznik. Berlin 1933. 

4. Taylor, G. I. : An Experiment on the Stability of Superposed Streams 

of Fluid. Proc. Cambr. Phil. Soc., vol. 23, p. 730, 1927* 

5. Schmidt, W. : Der Massenaustausch in freier Luft und vervandte 

Erscheinungen (Probleme der kosmischen Physik). Hamburg 1925* 

. 11 

6. Prandtl^ L. : Einfluss stabilisierender Krafte auf die Turbulenz. 

Vortrage ans dem Gebiete der Aerodynamik und verwandter Gebiete. 
Aachen 1929- p. 1. 

IT 

7. Prandtl, L. : Uber die ausgebildete Turbulenz. Vhdlg. des II. 

Intern. Kongr. f. techn. Mechanik, Zurich 1926. 

8. Prandtl, L. : Meteorologlsche Anwendung der Streamings lehre. 

"Beitr'&ge zur Physik der freien Atmosphere" (Bjerknes— Festschrift) , 
vol. 19, p. 188 , 1932. 

9 . Richardson, L. F. : The Supply of Energy from and to Atmospheric 

Eddies. Proc. Roy. Soc. London, vol. 97; P- 354 (1920). 

10. Taylor, G. I. : Internal Waves and Turbulence in a Fluid of Variable 

Density. Rapp. Proc.— Verb. Cons. Intern, p. l’Expl. de la Mer, 
LXVTI. Kopenhagen, pp. 20 to 30 (1931)- 

11. Taylor, G. I. : Effect of Variation in Density on the Stability of 

Superposed Streams of Fluid. Proc. Roy. Soc. A. Vol. 132, p. 49 9; 
1931. 

12. Goldstein, S. : On the Stability of Superposed Streams of Fluids of 

Different Densities. Proc. Roy. Soc. A. Vol. 132, p. 524, 1931. 

13. Tollmien, W. : Uber die Entstehung der Turbulenz. Nachr. d. Ges. d. 

Wiss. zu Gottingen, Math.— Phys. Klasse 1929; P- 21, and Verhandlg. d. 

III. Intern. Kongr. f. techn. Mech. Stockhom 1930, p. 105. 



NACA TM 1262 


47 


14. Rayleigh, Lord: On the Stability or Instability of Certain Fluid 

Motions III. Proc. London Math. Soc. 27, 18 95. 

15. Sommerfeld, A.: ^Ein Beitrag jsut hydrodynamischen Erklairung der 

turbulenten Flussigkeitsbevegung. Atti d. IV. congr. int. del Mathem. 
Rome 1909. 

16. Mises, R. v. : Beitrag zum Oszillationstheorem, Heinrich Weber- 

Festschrift 1912$ derselbe: Zur Turbulenztheorie, Jahresber. d. 

deutsch. Mathem. Ver. 1912. 

17. Hopf, L. : Der Verlauf kleiner Schwingungen in einer Streaming 

reibender Flussigkeit. Ann. d. Phys. 44, p. 1, 1914, and: Zur 
Theorie der Turbulenz. Ann. d. Phys. 59# p. 538, 1919* 

18 . Prandtl, L. : Bemerkungen Uber die Entstehung der Turbulenz. Ztschr. 

f. angev. Math. u. Mech. vol. 1, p. 431, 1921. 

19. Tietjens, 0.: Beitr^ge zum Turbulenzproblem. Dias. Gottingen 1922 

and. Ztschr. f. angew. Math. u. Mech. vol. 5# P. 200, 1915- 

20. Wendt, F. : Turbulent e Stromungen zwischen zwei rotierenden koaxialen 

Zylindem. Ing.— Arch. vol. IV, p. 577# 1933- 

21. Schlichting, H. : ffber die Entstehung der Turbulenz in einem 

rotierenden Zylinder. Nachr. d. Ges. d. Wiss. zu Gd&t ingen, Math.— 
Phys. Klasse 1932, p. l60. 

22. Schlichting, H. : Zur Entstehung der Turbulenz bei der Plattenstromung. 

Nachr. d. Ges. d. Wiss. zu Gottingen, Math.-Phys. Klasse, 1933# 
p. l8l. also: Ztschr. f. angew. Math. u. Mech. vol. 13# 1933- 

23 . Blasius, H. : Grenzschichten in Flussigkeiten mit kleiner Reibung. 

Diss. Gottingen 1907 and Ztschr. f. Math. u. Phys. vol. 51# p. 1# 

1908. 

24. Prandtl, L. : Uber Flussigkeitsbewegung bei sehr kleiner Reibung. 

Verhdlg. d. III. Int. Math.-Kongr. Heidelberg 1904, or "Vier 
Abhandlungen zur Hydrodynamik und Aerodynamik, Gottingen 1927. 

25. Prandtl, L. : Modellversuche und theoretische Studien uber die 

Turbulenz einer geschichteten Luftstronning. Deutsche Forschung 
No. 14, Berlin 1930. 

26. Prandtl, L. , and Reichardt, H. : Einfluss von Warmes chi chtung auf 

die Eigenachaften einer turbulenten StrSmung. Deutsche Forschung 
No. 21, Berlin 1934. 


48 


MCA TM 1262 


B/asius profile 
Approximation function 


y/t 


1.0 


1 

i 

i 

>1 

i 

i 


\ 

Temperature A 
\ 

Velocity — 


k / 


yr 

\ 

y' 5 

\ 

0 T °C \oo 


0.5 U/Um t.O 


Figure 1.- The investigated laminar flow with stratification of densities. 



Figure 2.- Polar diagram 
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Figure 7.- The critical Reynolds number of the plate flow as a function of 


the Richardson number and the Froude number = R^ ( ©, 
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Figure 8.- Stability limit 0^ as a function of Froude number — 
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Figure 14. 



Figure 15.- Comparison with test data. In agreement with theory, oi^y 
laminar states were observed in the stable range, and only turbulent 
states in the unstable one. 
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